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The generalized Master equation of the Nakajima-Zwanzig (NZ) type has been used extensively 
to investigate the coherence dynamics of the central spin model with the nuclear bath in a narrowed 
state characterized by a well defined value of the Overhauser field. We revisit the perturbative NZ 
approach and apply it to the exactly solvable case of a system with uniform hyperfine couplings. This 
is motivated by the fact that the effective Hamiltonian-based theory suggests that the dynamics of 
the realistic system at low magnetic fields and short times can be mapped onto the uniform coupling 
model. We show that the standard NZ approach fails to reproduce the exact solution of this model 
beyond very short times, while the effective Hamiltonian calculation agrees very well with the exact 
result on timescales during which most of the coherence is lost. Our key finding is that in order 
to extend the timescale of applicability of the NZ approach in this case, instead of using a single 
projection operator one has to use a set of correlated projection operators which properly reflect the 
symmetries of the problem and greatly improve the convergence of the theory. This suggests that the 
correlated projection operators are crucial for a proper description of narrowed state free induction 
decay at short times and low magnetic fields. Our results thus provide important insights toward the 
development of a more complete theory of central spin decoherence applicable in a broader regime 
of timescales and magnetic fields. 



I. INTRODUCTION 

Qubits based on spins of single electrons confined in 
a semiconductor environment* 1 ^ are one of most promis- 
ing platforms for quantum computation. When the host 
semiconductor is a III-V compound material (i.e. GaAs 
or InAs), the presence of the nuclear spins, coupled to 
the electron (or a hole) by hyperfine (hf) interaction, is 
unavoidable £1 The qubit decoherence induced by hf in- 
teraction with such a nuclear bath has been a subject 
of many theoretical works. For this paper, the directly 
relevant ones are those focusing on purely hf-induced 
dynamics^ - — i.e. on the problem described by the cen- 
tral spin Hamiltonian 

h = Y / A kS-i k + ns z + Y,ukik , (i) 

with S denoting the central spin, the nuclear spins, 
Ak the hf couplings, and fl and ujk the central spin and 
nuclear spin Zeeman energies, respectively. This Hamil- 
tonian can be solved exactly via Bcthe ansata±2i22i22 for 
the case of uniform ujk = cj, where U) can be removed 
by going to a rotating frame. 9 However, extracting po- 
tentially experimentally relevant information about the 
dynamics of the system, for example the time depen- 
dence of the reduced density matrix of the electron spin 
after initializing it in a pure state, is very hard when 
starting from the exact eigenstates ; 19 ' 30 ' 34 and the largest 
system sizes considered had up to only N = 30 nuclear 
spins. 19 Numerical calculations of the system's evolution 
can be also obtained using the Chebyshev polynomial 
method^i^ which has been used to calculate the free 



evolution^ and spin echo2£ signals for N <20. Larger 
systems with A^IO 4 were investigated using an approx- 
imate spin-coherent-state P representation , 14 i 35 which, 
however, can only be used in the case of an unpolar- 
ized thermal nuclear bath, which is too restrictive for 
many applications including the considerations presented 
here. In order to calculate the dynamics of the cen- 
tral spin on experimentally relevant timescales, and for 
experimentally relevant system sizes (N ~ 10 4 — 10 6 
in III-V compound quantum dots), one has to resort 
to approximate analytical methods, which include the 
use of the Generalized Master Equation (GME) of the 
Nakajima-Zwanzig (NZ) type^^^ the GME of the 
time-convolutionless (TCL) kind ) 22 ' 24 ' 36 and the cluster- 
expansion type theories using the effective pure de- 
phasing Hamiltonian ;^ ' 12 i 25 ~ 27 i 36 obtained perturbatively 
from the original Hamiltonian of Eq. ([T]) by an appropri- 
ate canonical transformation. 

In this paper we focus on the NZ theory applied to 
the case of narrowed state free induction decay (NFID), 
i.e. on the situation in which the nuclear system is 
prepared in an eigenstate of the Overhauser operator 
h z = J^k^-klk- This case is both important for quan- 
tum computation applications and experimentally rele- 
vant, since significant progress in nuclear state narrow- 
ing has been made recently 3 ^— (see also Ref. d and ref- 
erences therein). Furthermore, the high magnetic field 
regime of NFID (defined by the condition of fl > A where 
A = Y^,k A k ~ 90 fieV in GaAs) was extensively inves- 
tigated using the NZ approach i 9 ' 23 i 28 Let us note that 
the other GME approaches, such as the TCL formal- 
ism and, more importantly in the context of this paper, 
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the master equations employing the correlated projection 
operators ; 22 ' 24 ' 36 have been until now only employed in 
the case of thermal (non-narrowed) nuclear baths. 

Let us briefly recount the main results of the NZ the- 
ory for the case of NFID. In the NZ approach, the equa- 
tion of motion for the reduced density matrix of the cen- 
tral spin has an integro-differential for m 42 ' 43 with the 
memory kernel explicitly showing the influence of the 
history of the spin on its dynamics at a given point in 
time (i.e. the non-Markovian effects) . While this feature 
is physically appealing and offers a natural setting for 
investigation of the Markovian approximation and the 
corrections to it arising from the temporal nonlocality 
of the memory kernel, the NZ approach is technically 
very demanding. In the case of the central spin model, 
the memory kernel is expanded in powers of the flip- 
flop operator, Vg ~ S ± I T , and the expansion becomes 
very complicated with increasing order of perturbation 
theory. It also does not have any diagrammatic struc- 
ture since it is not based on expanding a generalized 
(i.e. time- or contour-ordered) exponential, which would 
allow one to write down more easily higher-order contri- 
butions and identify the most important classes of terms 
at each order. Furthermore, the integro-differential equa- 
tion is solved via Laplace transform, and the final result 
is obtained by inverse transform of a rather complicated 
expression. 

Despite these obstacles, large progress has been made 
with this method. The 2nd order solution- was shown to 
reproduce the exact solution for the fully polarized bath 
case* 7 - signifying the important role played by nuclear po- 
larization (which limits the phase space for flip-flops) in 
improving the convergence of the NZ approach. For a 
general polarization, the 2nd order result was shown to 
lead to a very small coherence decay at high B fields sat- 
isfying _4/i7<C 1. Specifically, the fraction of coherence 
lost was only about A 2 /NO 2 <C 1. A complete decay of 
the transverse electron spin components was obtained af- 
ter going to the 4th order of the flip-flop expansioni 23 i 28 
There, a controlled solution was obtained for A/0< I in 
the case of an unpolarized bath: after an initial quadratic 
decay shoulder the decay was of the exponential form 
exp(-i/T 2 ), with T 2 - NO 2 /A 3 , with the l/t 2 tail ap- 
pearing at very long times. Recently it was showr* 2 ^ 
that with finite nuclear polarization, and for nuclear spin 
/ < 1, it is possible to extend the 4th order solution to 
lower B fields (down to A/O, ~ 1 for unpolarized nuclei). 
As O is decreased, corrections to the above formula for 
T 2 appear, and at A/ Ok 1 a minimum of T 2 and a new 
kind of coherence envelope oscillation were predicted^ 

It has to be stressed that most of the above fea- 
tures of the high-field NZ solution, especially the ex- 
ponential decay, arise due to the inhomogeneous cou- 
pling of the nuclei to the central spin. In the 4th or- 
der in electron-nuclear flip-flops, the processes in which 
pairs of remote nuclei flip-flop among each other con- 
tribute to the system's dynamics. In the effective Hamil- 
tonian language ; 12 ' 16 ' 23 ' 25 ' 26 one can say that the Over- 



hauser field fluctuates due to electron-mediated flip- 
flops of these remote nuclei. At short times defined 
by t < N/A — 1/A k the differences of Knight shifts 
Aki = Ak — Ai of these nuclei are unimportant due to 
the energy-time uncertainty, leading to the independence 
of the quadratic decay shoulder on the wavefunction 
shaped At longer times, the Knight shifts must ful- 
fill more stringent energy conservation conditions, giv- 
ing rise to the long-time Markovian exponential decay, 
with the T 2 time strongly dependent on the shape of 
the wavefunctioni^ 3 - (Compare also the expressions for T 2 
from Refs. [H and [26| which correspond, respectively, to 
the two-dimensional Gaussian e~ [x +y >' a and the same 
function modified by a cosine form factor in the z direc- 
tion). 

An insight into the importance of the short-time 
regime is provided by theories based on an effective- 
Hamiltonian, specifically the Ring Diagram Theory 
(RDT), in which a class of diagrams of leading order 
in 1 /N is resummed in the cluster expansion of the elec- 
tron's density matrixi 4 - ' 25 ' 26 According to these theories, 
for low fields such that A/VN O < A, most of the 
coherence decay occurs at short times, at which the dis- 
tribution of the hf couplings does not matter, and only 
the quantities A and N are important. This suggests 
that the low-field dynamics of NFID in a realistic system 
is closely related to dynamics in a system with uniform 
hf couplings. 

In RDT, one starts with an approximate effective 
Hamiltonian derived in the 2nd order of perturbation the- 
ory with respect to the flip-flop operator! 11 ' 12 ' 23 ' 26 While 
this step makes it harder to precisely ascertain the lim- 
its of applicability of this approach, it allows for a very 
convenient formal simplification of the problem. The ef- 
fective Hamiltonian is of the pure dephasing form, which 
allows for the use of standard tools of non-equilibrium 
quantum dynamics, such as the closed time-loop contour 
and the cluster expansion of the bath average of a gen- 
eralized exponent i 17 ' 25 ' 26 Such methods were previously 
used in order to calculate decoherence for many other 
models of quantum baths,— ~— however in the spin bath 
case there are additional complications due to the lack of 
a simple Wick's theorem for spin operators. (See Refs.[l7l 
and |47| for the derivation of all the diagrams in the 4th 
order of perturbation theory with respect to inter- nuclear 
interactions.) However, for the electron-mediated (or hf- 
mediated) interactions, due to the fact that they couple 
all the N nuclear spins with comparable strength, one 
can resum all the leading diagrams in the 1/N expan- 
sion, which amounts to calculating only one diagram at 
each order of the cluster expansion. 

For NFID at large B fields (A/0<^1), this resumma- 
tion leads to an exponential decay with the same T 2 as the 
one obtained from the NZ approach 23,28 (albeit without 
the the l/t 2 long-time tail), down to A/Om 1 at which 
the minimum of T 2 is not reproduced. Within RDT one 
can calculate the short-time decoherence also at low mag- 
netic fields obeying only the condition A/VNO <C 1. 
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The spin echo signal predicted by RD T 25 i 26 was later 
seen in experiments on double GaAs quantum dots, 48 
and it was recently rederived using an explicitly semi- 
classical approach^ For NFID, a decay envelope of the 
form (l + ^/r) 2 )- 1 / 2 , with t^NQ/A 2 , and a tt/2 phase 
shift of the electron precession for t^r were predicted^ 
As we show explicitly in this paper, this result arises also 
as the large N limit of the exact solution of the fully 
quantum calculation for the system with uniform hf cou- 
plings. 

In this paper we apply the standard 4th order NZ 
theor y 9 ' 23 ! 28 to the otherwise exactly solvable model with 
uniform hf couplings. We find that the NZ calculation 
of NFID disagrees with the exact result beyond a very 
short timescale (at which the coherence barely decays), 
while the RDT calculation reproduces the envelope of the 
NFID signal very well. We show that in order to improve 
the performance of the NZ theory, one needs to replace 
the commonly used single projection operator by a fam- 
ily of projection operators on nuclear subspaces which 
are singled out by the electron-nuclear coupling. Such 
an approach was used previously in both the NZ and 
TCL generalized master equation theories in the simpler 
case of a thermal nuclear bath; 22 i 24 ' 36 but it has not been 
used until now for the narrowed FID case. This result 
strongly suggests that the main physics of the short-time 
and long-time central spin decoherence is significantly 
different, with simple dephasing by distinct nuclear states 
dominating at short times, and nontrivial dynamics of in- 
homogeneously coupled nuclei dominating at long times. 
The use of the correlated projection operators allows the 
NZ theory to fully capture the short-time dynamics in a 
wide range of magnetic fields. 

Let us briefly note that there exists a family of theo- 
ries investigating the problem of dipolarly-induced elec- 
tron spin dephasing , 12 i 17 ' 47 i 49 ~ — in which the dipolarly- 
induced flip-flops between the nuclear spins lead to fluc- 
tuations of the Overhauser field. These theories, which 
are all based on some version of a linked cluster (cumu- 
lant expansion) theorem, are in agreement with the spin 
echo measurements in phosphorus doped silicon, 54 bis- 
muth doped silicon, 55 and in GaAs singlet-triplet qubits 
at magnetic fields higher than B 0.3 T4£ The dipolar 
processes are indeed expected to dominate the decoher- 
ence at high B since the increasing qubit energy splitting 
is suppressing the dynamical processes due to hf interac- 
tion, while the dipolar flip-flops of the nuclei are practi- 
cally unaffected by B. These experiments suggest that 
the window of parameters in which pure hf interactions 
dominate the electron spin decoherence is rather small 
in GaAs and Si. In fact, the spin echo results in large 
GaAs dots^ suggest that purely hf effects could be seen 
in lateral GaAs dots only at B<0.3 T. However, it was 
suggested that in small and strained InGaAs quantum 
dots the dipolar interactions between the nuclei might be 
suppressed by strongly inhomogeneous Knight shifts and 
by the quadrupolar interactions. 28 Experimental results 
suggesting suppression of dipolarly-induced spin diffusion 



in such dots have appeared recently^ 

Nothwithstanding the future experimental develop- 
ments, which hopefully will allow for measurements of 
purely hf-induced decoherence in a wide range of B fields, 
spanning the short-time semiclassical and long-time 
Markovian regimes (with possible nontrivial crossover be- 
tween them suggested by results of Ref . [H) , the central 
spin decoherence problem remains theoretically interest- 
ing in itself. Its further investigation will hopefully lead 
to establishing deeper connections between various ap- 
proaches known from the theory of open quantum sys- 
tems, and the properties of the exact Bethe ansatz solu- 
tion. 

The paper is organized in the following way. In Section 
HTl we introduce the hyperfine Hamiltonian and give the 
explicit form of the initial nuclear density matrix charac- 
terizing NFID. In Section UTTl we present the exact solu- 
tion of NFID for the uniform coupling model. Section ITVl 
offers an extensive review of the Nakajima-Zwanzig Mas- 
ter equation treatment of the central spin model, with 
the electron spin coherence computed up to 4th order in 
the hyperfine flip-flop interaction. This section, together 
with detailed appendices |A1 and IB1 can be mostly skipped 
by a reader deeply familiar with previous derivations of 
this theoryi^i 2 ^ However, we hope that our derivations, 
apart from making the paper self-contained, are useful 
additions to the existing works, clarifying some technical 
issues and expanding the discussion of certain aspects of 
a highly technical theory. In the subsequent Section fVl 
we present a derivation of NFID in the uniform coupling 
model using the 4th order NZ approach and show that 
the results disagree with the exact solution except at very 
early times. We demonstrate that the strong discrepancy 
between the exact and NZ solutions can be lifted by re- 
placing the standard single projection operator in the NZ 
theory with a family of correlated projection (CP) oper- 
ators, and that the modified NZ theory (which we refer 
to as the NZ-CP theory) then gives results which are in- 
distinguishable from the exact solution. We further show 
that the NZ-CP result is closely related to and improves 
upon that which is obtained by applying RDT to the 
uniform coupling model. The reasons for this behavior 
of the NZ theory, together with the implications of these 
findings for the coherence decay in realistic systems with 
non-uniform couplings, are discussed in Section |VI[ 



II. HAMILTONIAN AND INITIAL NUCLEAR 
BATH STATE 



The Hamiltonian for a single electron spin coupled to 
a nuclear spin bath via the hyperfine interaction has the 
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form 



H = Ho + Vg, 



with 



H a = ns z + J2^i z k + s z h z , 



v s = -(h+s- + h-s + ), 



(2) 
(3) 

(4) 
(5) 



The S l are the components of the electron spin operator, 
while the P k are the components of the spin operator for 
the fcth nucleus. The electron spin raising and lowering 
operators and their nuclear counterparts are defined 
in the usual way: 



S ± = S x ±iS y , 



It 



(G) 



We have included a magnetic field which points in the 
z-direction, and Q and ojk are the Zeeman energies of 
the electron and nuclei. In order to properly describe 
the nuclear bath in III-V semiconductors, we should al- 
low for several different isotopic species of nuclei. Rather 
than introduce an extra index to label distinct nuclear 
species, we have absorbed this index into the site index 
k on the hyperfine couplings Ak and nuclear Zeeman en- 
ergies uik ■ These quantities will in general depend on the 
nuclear species since different species will have different 
gyromagnetic factors. For example, the hyperfine cou- 
plings are given by Ak — voAa^^k)^ , where v is the 
volume of the primitive unit cell, ^(r^) is the electron en- 
velope wavefunction evaluated at the position of the 
fcth nucleus, and A a is the total hyperfine interaction en- 
ergy between the electron and a nucleus of species a and 
depends on the electron and nuclear spin gyromagnetic 
factors. 

The sums in Eqs. ©-(O range over some number N 
of nuclei which comprise the bath. Quantum dots re- 
siding in III-V materials like GaAs and InAs typically 
contain 10 4 -10 6 nuclei which interact appreciably with 
the electron, so physical values of N should lie in this 
range, although in analytical calculations often one sim- 
ply performs a sum over an infinite number of spins (out 
of which only about N are appreciably coupled to the 
electron spin). It is often useful to define the "effective 
number of nuclei appreciably interacting with the elec- 
tron" N: 



N 



A 2 



12k A k 



where 



A = J2^ k 



(7) 



(8) 



is the total hyperfine interaction energy between the elec- 
tron and the nuclear bath. For the most part, we will 



leave N arbitrary in subsequent sections, but when a 
more precise definition is necessary, we will use Eq. (|7J|. 

Let us quote some numbers in order to make a connec- 
tion to experimentally relevant timescales and magnetic 
field regimes. In GaAs, A~90 peV (see e.g. Ref. y| and 
references therein). Using the effective g- factor g e g^0.5 
of an electron in this material, Q « A for a magnetic field 
of B m 3 T. Another characteristic magnetic field in the 
central spin problem is that corresponding to the typi- 
cal Overhauser field for a thermal nuclear bath, A/V~N, 
which is on the order of a few mT. Note that the maximal 
hf coupling of a nuclear spin is Ak ~A/N, which is much 
smaller than Q for the magnetic fields considered in this 
paper (i.e. fields such that f2 > A/yN). An important 
quantity is also the timescale of N/A, defining here the 
boundary of the short-time regime. For N = 10 4 - I0 6 
this corresponds to timescales of 100 ns - 10 /us. 

It is difficult to obtain explicit results from the NZ 
GME unless strong assumptions are made about the form 
of the initial nuclear density matrix, which we denote 
pi(0). Fortunately, the NZ GME seems most tractable 
in the case of NFID, which is the primary focus of the 
present work. As we develop the structure of the NZ 
GME adapted to the central spin problem in Section HVl 
we will attempt to keep the form of pi(0) as general as 
possible at each stage of the calculation and make clear 
the points at which it is necessary to further specify pi(0). 
For the sake of clarity, we will state here the precise final 
form that we will impose on pi(0) in order to obtain ex- 
plicit results for the electron spin coherence in the NFID 
case. 

In discussing the initial density matrix of the nuclei, 
we have already alluded to the first assumption we will 
make about the initial conditions, namely that the initial 
density matrix for the total system is separable: 



p(0) =/> e (0)®p 7 (0). 



(9) 



Here, p e is the reduced density matrix describing the 
electron spin degrees of freedom. This assumption is 
quite reasonable since, in quantum computing applica- 
tions, it is generally desirable to initialize the electron 
spin to some well defined state which is independent of 
the nuclei. 

In contrast, a much stronger assumption we will make 
is that P/(0) is diagonal in the basis of h z eigenstates. 
As in Refs. HJH, we will denote these states by \rii) with 
i = 1...2 N ; they are just tensor products of eigenstates of 
the nuclear spin operators 1%: 

= l /fc ' m k) > J fe \lk, m D = m l \!k,ml) , 
fe 

(10) 

where Ik is the spin of the fcth nucleus. We may write 

h z \n i ) = hl i \n i ), h z n ^J2 A ^- ( n ) 

fe 

We will further assume that pi(0) only contains compo- 
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Pi{°) = ^2 Pa l n *> 



where 



nents which share the same h z eigenvalue, denoted h z n : to the problem. In particular, this limit is exactly 

solvable ; 7 ' 13 ' 20 ' 22 ' 27 ' 57 and we will begin by reviewing this 
solution for the NFID case. In Section [Vj we will use 
this solution to test how well the NZ approach (both a 
standard one, and the one using correlated projection 
operators) can describe the dynamics of the box model. 

Restricting attention to the homonuclear case, Uk = w, 
with all the nuclei having spin / = 1/2, we can solve the 
uniform coupling model exactly by working in the basis 
of eigenstates of the total nuclear spin operator I z = 
J2 k Ik, we denote these states by \jm) with I z \jm) = 
m\jm). In this basis, we may write the expectation value 
(S + (t)} in terms of the appropriate matrix element of the 
reduced density matrix evolved with respect to the full 
Hamiltonian: 



hC 



Vt e l...g, 



(12) 



(13) 



The form of the density matrix given in (|12M13[) describes 
a particular "narrowed" set of allowed states for the nu- 
clear bathi2 In the so-called "box model" limit where the 
hyperfme couplings are all equal, Ak = A = A/N, this 
amounts to restricting the possible nuclear states to a set 
which contains only states with a fixed net polarization 
along the magnetic field direction. We will see that the 
NZ GME approach enjoys many simplifications when the 
nuclear bath states are constrained in this way. 

For some of the more exp licit results we will obtain, we 
will also follow Refs. I9ll28l in assuming that the nuclear 
spin bath is uniformly polarized. By this, we mean that 
all traces of nuclear spin operators are independent of the 
nuclear site indices. For example, the correlator 



Tr{l+^ P/ (0)/ fe -/+} 



(14) 



is independent of k,£,p,q under this supposition. The 
uniform polarization assumption is less crucial than the 
narrowed-state condition of Eq. (|12|) and can be lifted 
in most cases, albeit at the expense of having to deal 
with more complicated algebraic expressions. For the 
purposes of elucidating more general features of the NZ 
GME approach, it is often a useful assumption to make. 
Throughout this paper we make explicit the points at 
which we invoke the uniform polarization condition. 



III. UNIFORM COUPLING MODEL 

In this section, we will consider the uniform-coupling 
(also referred to as 'box model') limit in which all the 
hyperfme couplings are equal, 



A k = A = 



A 

N 



(15) 



(S + (t)) = 



YV- (|,j,m|e- 4i/? - 



p(0y Ht \t,j,m). (16) 



The f, 4 inside the kets denote the eigenstates of the elec- 
tron spin operator S z . p(0) is the full density matrix at 
time t = 0, and tij is the degeneracy of nuclear states 
with a fixed j and m, 58 



2.7 + 1 



(f-i)Kf +i)!f + 



(17) 



In order to facilitate the computation, we will assume 
that the initial nuclear density matrix is diagonal in the 
\jm) basis. In fact, in order to compare with the NZ 
GME result computed in Sec. [V] we will assume that 
the initial nuclear density matrix is proportional to the 
identity in a subspace of fixed m and is zero outside this 
subspace. This is equivalent to the form of pi(0) that 
we assumed in Eqs. (|12[) -(fT"3 j) . if in the latter we make 
the further assumption that the nuclei are uniformly po- 
larized, pa = 1/Z with Z = Y^j n j- When the density 
matrix has this form, our formula becomes 



This is a semiclassical limit of the central spin problem 
since we have a large collective nuclear spin I = ^ fc 1^. 
coupled to the electron spin — the dynamics is mostly 
classical for large N. This is also very close to the static 
bath limit: for Q A/\/~N and for a typical collective 
spin magnitude |/| ~ V~N, I is almost static as it cannot 
follow the quickly precessing electron spin vector (if we 
neglect the nuclear Larmor precession, which is crucial 
for the case of spin echo decay 25 ' 26 ' 31 ' 48 but which is of 
much smaller importance for NFID). 

The limit of uniform hyperfme couplings greatly sim- 
plifies the electron-nuclear dynamics, making this a 
natural limit in which to compare various approaches 



N/2 



(S + (t)) = \ £ rij (i,j,m\e- iHt Pe(0)e im |t, j,m) 



where p e (0) is the initial electron density matrix. 



(18) 



In the case where all hyperfme couplings are equal, 
Ak = A, it is straightforward to work out the action of the 
evolution operator e~ lHt on the |f / \.,j,rri) states; 22 ' 27 



-iHt 



-iHt 



It, J, m) 
\i,j,m) 



|t,J,m) + b jm \l,j,m+ 1) , 



It, j, m -1) ,(19) 



6 



with 



iEZt 



COS 



N 



jm 



m 



X + ( N + 



-te 



-iEZt jm 

N+ 



jm 



iEZt 



COS 



N 



Jin 



iEZt jm 



N 



jm . 



N 



(20) 



jm 



and 



V ■ = 

jm 

Z m = 



7V± 



[m± 1/2)lu - A/(4N), 
Ay/j(j + 1) - m(m ± 1)/N, 
±[n-w + A(m±l/2)/N], 

yj(xf m y + (z±)*. 



(21) 



Using (|19p and expanding the initial electron density ma- 
trix as 

p e (0) = U+2(S Z (0))S Z + (S+(0))S- + (S-(0))S+, (22) 

one finds that only the S 1- component of p e (0) con- 
tributes to {S+(t)): 



(S + (0)) 

z 

<S+(0)) 



JV/2 

J2 n i a > i,m|e- <Ht S'- e ifft |t > i,'. 

j=|m| 
JV/2 



E 1 

3=M 



(23) 



For later comparison, we transform this result to a certain 
rotating frame defined with respect to the frequency Q n + 
AQ, where 



n„ = n + hi 



(24) 



and the "Lamb shift" Ail will be defined later on. In the 
rotating frame, the exact solution to the uniform coupling 
model is then given by 

N/2 

•Eexactit) — ^ / ^ ^ ^j^jm^jni- (^^) 

j=\m\ 

To obtain this expression, we have used that h z n = Am/N 
in the case of uniform coupling. 



parts, often referred to as the relevant and the irrelevant 
parts 



.42.43 



P = Prel + Pi 



(26) 



The relevant part, p re i, is a density matrix for the degrees 
of freedom whose evolution one wishes to compute, while 
the irrelevant part, pi rr , contains the remaining degrees 
of freedom in the total system. While these remaining 
degrees of freedom will influence the evolution of p re i, 
their own evolution is not of direct concern and need 
not be computed explicitly. In applications where one 
considers the dynamics of a system coupled to a bath, 
p re i is typically a density matrix describing the degrees 
of freedom of the system. In this case, note that strictly 
speaking, pi rr is not a density matrix since it does not 
satisfy Trp^ rr = 1. 

The partition described in Eq. (|26[) is implemented by 
introducing projection superoperators P and Q such that 

Pp = prel, Qp = Pirr, P + Q = 1, PQ = 0. (27) 

The last two relations above ensure that P and Q are 
true projectors in the sense that P 2 — P and Q 2 = Q. 
With these operators, the Liouville equation for p can be 
transformed into an exact equation for the evolution of 
p re i , which in the case of a time- independent Hamiltonian 
H has the form^d3 

Pp{t) = -iPLPp(t) - i [ dt't(t - t')Pp(t'), (28) 



T,(t) = -iPLQe~ tLQt QLP. 



(29) 



The Liouvillian superoperator L implements the evolu- 
tion of the total system and is defined to act on an arbi- 
trary operator O according to 



LO = [H,0\. 



(30) 



The superoperator S is referred to as the memory kernel; 
this quantity contains the full dynamics of the bath and 
controls how these dynamics affect the evolution of the 
system. It is important to note that the particular form 
of the NZ equation given in Eq. (|28|) assumes that the 
initial density matrix satisfies the condition Qp(0) = 0. 

In the case of our central spin model, Eqs. ©-([S]), we 
choose prei to be essentially the reduced density matrix 
for the electron spin, p e . This means that the projection 
operator P involves a trace over the nuclear bath£ i 23 ' 28 



Pp = P I (0)Ty iP = P e® pi{Q). 



(31) 



IV. NAKAJIMA-ZWANZIG MASTER 
EQUATION AND ITS PERTURB ATIVE 
EXPANSION 

The first step in constructing the NZ GME is to sep- 
arate the density matrix for the total system into two 



We have introduced a factor of the initial nuclear density 
matrix pi(0) in the definition of P to satisfy the con- 
traint P 2 = P. Note that the condition Qp(0) = is 
tantamount to assuming that p(0) is separable: 



p(0) =p e (0)®p/(0). 



(32) 
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This is an assumption we will make throughout this work. 

It is important to stress that Eq. (f3~Tj) is not the only 
possible choice for the projector P. For example, it is 
possible to instead define P as a sum over many projec- 
tion operators which project onto various subspaces of 
the nuclear bath state space, and the choice of P can 
strongly influence the convergence properties of the re- 
sulting theoryi 22 i 24 In fact, we will see later on in the 
context of the uniform coupling model that the choice 
made in Eq. pip is far from ideal, and that a more 
sophisticated choice leads to a vast improvement in the 
convergence of the theory. Given the connection between 
the box model and the short-time, low B-field regime of 
the real system with non-uniform couplings, we will ar- 
gue that this observation carries important consequences 
for the development of a more complete theory of the 
central spin problem. For now, we will keep the defini- 
tion of P in Eq. pip and continue our construction of 
the "standard" NZ theory. 

The form of the NZ equation given in Eq. (f^5|) is quite 
difficult to work with. In order to reduce it to a more 
tractable form, we will place restrictions on the structure 
of the initial nuclear density matrix, pi(0). In particular, 
we will assume that pi(0) is diagonal in the |nj) basis: 
pi(0) = pu \n.i) {fii\- If we then multiply both sides 
of Eq. (|28p by the operator S + and take the trace, we 
obtain an equation for (S + ) only: 



d_ 

di' 



(5+(t)> =iSl n (S + (t))-i / dt"E(t-t')(S + (t')), (33) 



where O n = Q + hZ with 



P,i 



(34) 



Note that this definition of h z n is more general than the 
meaning we gave this symbol in Eq. (I13p in the context 
of the fully restricted nuclear density matrix. The two 
meanings coincide when pi(0) has the form of Eq. (fT12j). 
The memory kernel is now a function instead of an oper- 
ator: 

E(t) = -iTr [S + PLQe~ iLQt QLPS- Pl {0)] . (35) 

Without the assumption that pi(0) is diagonal in the \rii) 
basis, additional terms involving (S z (t)) and (S~(t)) ap- 
pear in Eq. (|33|) , and the problem becomes considerably 
more complicated. Note that we did not have to assume 
the fully restricted form for pi(0) quoted in Eq. ([12")) 
where the sum is only over \m) states corresponding to 
the same Overhauser field. Eq. (|33|) is valid when the 
sum contains other states as well, i.e. when the nuclear 
state is not narrowed. 

Eq. p"3"|) is an integro-differential equation which can 
easily be solved by performing a Laplace transform, after 
which the equation becomes algebraic with the solution 



dte- st (S + (t)) 



(s + (t = o)} 

— i£l n + iS(s) 



(36) 



The solution in the time domain is then obtained by com- 
puting the Bromwich inversion integral, 



(S + (t)) = 



1 

2iri 



7+ioo 



dse st (S + {s)), 



(37) 



where the contour defined by the real number 7 must be 
chosen such that it lies to the right of all the poles of 
(S + (s)). Therefore, solving for (S + (t)) requires solving 
for the Laplace transform of the memory kernel: 



S(s) = / dte~ st E{t) 
Jo 



-iTr 



S + PLQ 



1 



s + iLQ 



QLPS-pi{0) 



(38) 



Computing the memory kernel exactly is a difficult prob- 
lem except in very simple cases, so we will proceed to 
calculate it perturbatively in the following section. 



A. The expansion of the memory kernel 

In order to make further progress, it is necessary to 
construct a perturbative expansion of the memory kernel 
given in Eq. (1381) . We will expand this quantity in a 
power series in the number of flip-flops, i.e. in powers of 



S( S ) = E( 2 »( S ) + E( 4 >( S ) + 0(y ff 6 ) 



(39) 



We can only have even terms in this expansion since the 
flip-flops are virtual in the physical limit Q 3> u>k- For 
nuclear density matrices pi(0) which are diagonal in the 
\rii) basis, the odd terms are strictly zero due to the struc- 
ture of the memory kernel, Eq. (|38|) . We will see shortly 
why a zeroth order term has not been included in Eq. 
(EH). 

To facilitate the expansion, it is convenient to define 
the following superoperators: 



L O = [H ,O], 



L v O = [V s ,0\. 



(40) 



(41) 



Replacing L — >• Lq + Ly in (|38|) . we can then recast our 
perturbative flip-flop expansion as an expansion in Ly. 
The expansion of the memory kernel is tantamount to an 
expansion of the operator s+ ^ Lo 1 fLi/ ^Q in powers of Ly, 
and this series is straightforward to construct: 



{1 - iG Q {s)L v Q ~ Gq(s)L v QGq(s)L v Q 

(42) 



s + iLQ 

+iG Q {s)L v QG Q {s)L v QG Q {s)L v Q}G Q {s) + 0(L V ), 



where 



Gq(s) 



iLo( 



(43) 
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Using this expansion, we will compute the memory kernel 
up to fourth order. 

Before we proceed to compute the terms of the series, 
we first pause to write down some identities which will 
be useful at various stages of the expansion. When pi(0) 
is diagonal in the |n.j) basis, the following identity holds 
for any integers k and £: 



PL«L V L* P 



0. 



(44) 



This identity can be further generalized to include any 
combination of L 's and an odd number of Ly's sand- 
wiched between two P's. When p/(0) is further restricted 
to the form given in Eq. (|12j) , we have an additional iden- 
tity (expressed in two equivalent ways): 



QL P = 0, QL Q = QL Q . 



(45) 



The first form of this identity allows us to replace one of 
the L's in Eq. (J3HJ) with Ly. 



£(*) = -iTr[S + PLQ 



1 



s + iLQ 



QL V PS- Pl (0)i (46) 



which indicates that the zeroth order memory kernel van- 
ishes, as was already presumed in (|39j) . Note that in using 
Eq. (|45| to arrive at Eq. (j46|) . we are assuming the fully 
restricted form of p/(0) stated in Eq. (fT2|) -([T5 1) . whereas 
up until now, we have only needed to assume that pi(0) 
is diagonal in the |n;) basis. In fact, at this point it is not 
really necessary to impose Eq. (fT2|) , and we could instead 
continue to suppose only that pi(0) is diagonal, in which 
case we would find a non- vanishing zeroth-order term in 
the memory kernel expansion. However, slightly further 
into the calculation we will impose Eqs. (fl"2j) - (fT3"|) . and 
so at this stage we may as well make the simplifications 
that this form of pi(0) brings. 



B. Second-order memory kernel 

Inserting Eq. (l4"2l into Eq. (HU), we find at second 
order 

£( 2 >( s ) = -i r Tr[S + PL v QG Q (s)QL v PS-p I (0)] 

Tr[5+ ■PL QGq(s)L v QGq(s)QL v PS- Pi (0)]. 

(47) 

The second term on the right-hand side can be shown to 
vanish identically for the central spin Hamiltonian. Fo- 
cusing then on the first term, this expression is somewhat 
complicated by the dependence of Gq(s) on Q. Instead 
of working with Gq(s), we choose to work with its time- 
domain counterpart, e ~ lL oQ t j m terms of which we have 

E (2) (i) = -iTr[S + PL v Qe- lLoQt QL v PS- pj(0)]. (48) 

It is easy to show that the identity ([4"5"j) enables the ex- 
traction of the projection operator Q from the exponent: 



-iLoQt 



-iLot 



QG(t). 



(49) 



The superoperator G(t) is simply the evolution operator 
corresponding to the unperturbed part of the Hamilto- 
nian, Hq\ 



G{t)p(0) = e- l " ot p(0)e l 
Also notice that (|49[) immediately implies 

QG Q (s) = QG(s), 
with the "propagator" defined as 

1 



G(s) = 



iLn 



(50) 



(51) 



(52) 



At this point, we will invoke (by using Eq. (I45[l ) the most 
constrained form of the initial nuclear density matrix, 
Eq. p^)) - ([T5| . as it is not clear how to proceed without 
factoring the projector Q out of the propagator G(s). If 
one could proceed without performing this factorization, 
then it would suffice to assume only that pi(0) is di- 
agonal in the \rii) basis, albeit several more terms would 
contribute at each order of the memory kernel expansion. 
In any case, we are ultimately interested in applying this 
formalism to the case of NFID, which is defined precisely 
by Eqs. dH-gSJ. 

Plugging Eq. (|4T))) into Eq. (|4"5)l . eliminating factors 
of Q with applications of the identity given in Eq. (|44l) , 
and dropping factors of P using the observations that 
PS-pi(0) = 5-/9/(0) and Tr[S+PO] = Tr[5+C] for any 
operator O, we arrive at 

£ (2 >(t) = -iTr[S + L v G(t)L v S~ pi(0)}. (53) 

It remains to plug in explicit expressions for the various 
operators and perform the trace. Defining the following 
set of operators which act in the nuclear subspace, 

U±(t) = e ^i( n +^)*- i Sfc^ J fc* ) (54) 

we can express the action of G(t) on a matrix p as 



Go = ( U +Pn u + u +Pu ui 
\u-pxul U- Pu u!_ 



while the action of Ly has the explicit form 



L v p 



1 fh~ 

2 \h+p n 



Pit 



Ptl h ^ 
Pll h ^ 



h Pll 

h + Ptl 



p^h 
Plt h " 



(55) 



(56) 



where for instance pf± = (t\p\l)i with It) and \\) de- 
noting the eigenstates of the electron spin operator S z . 
After a bit of algebra, we find 



2 iwfct 



4/ 



cl e 



where 

4 = Tr{^/± P/ (0)} 
The Laplace transform of this is 



(57) 



(58) 



1 



S(2) (*) = ^E^- 



i(w fc - ^) s - i(uj k + 4»-) 



(59) 
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C. Lamb shift and rotating frame 

The effect of £< 2 ) on (S+(t)) was studied extensively 
in Ref. [§, where it was shown that it leads to both a shift 
in the precession frequency of the electron spin (Lamb 
shift) as well as to a small decay of (S + (t)) at short times 
referred to as "visibility loss" . In the limit of large N, the 
visibility loss effect is suppressed and one can make 
the approximation that the only role of the second-order 
memory kernel is to generate the Lamb shift. We will 
refer to this as the "Lamb shift approximation." 

To make the definition of the Lamb shift more precise, 
first recall the Bromwich integral formula from Eq. ([37| : 



(S + (t)) = 



(S+(t = 0)> 



2m 



7+ioo 



dse s 



1 



"f — ZOO 



- i£l n 4- z£(s) 

(60) 

Since E(s) vanishes at zero hyperfine coupling, it is 
clear from this formula that the zeroth-order behavior of 
(S + (t)) is just a precession with frequency il n . At non- 
zero coupling, this precession is shifted by an amount 
determined by the real part of the memory kernel. We 
define the Lamb shift, Ail, as the shift in precession fre- 
quency that would occur if the imaginary part of the 
memory kernel (the part which produces decay) were 
zero. We can read off a self-consistent equation for Ail 
from Eq. (|60|l by supposing Im[£(s)] = and requiring 
that the integrand have a pole at s = i(il n + Ail): 



Ail = -Re[£(ifi„ + iAfi)]. 



(61) 



In our flip-flop expansion, this can be approximated by 
Ail fa -Rc[£ (2) (^„ + iAil)]. (62) 

Using this formula in conjuction with Eq. (|59[) , we then 
find that the Lamb shift is given by 



(63) 



in the limit il n 3> ojk, Ak- 

Now that we have defined the Lamb shift more pre- 
cisely, we should also clarify the manner in which the 
Lamb shift approximation is implemented. This approx- 
imation consists of replacing Y,( 2 > (s+iil n +iAil) — > —Ail 
in the expansion of the memory kernel. In order to make 
this replacement in (|60p , we need to shift the integration 
variable s — > s + iil n + iAil. In the Lamb shift approxi- 
mation, we then have 



i(n„+AQ)i 

2iri 



7+zoo 



dse s 



1 



7 - loo s + i£( 4 )(s) 

(64) 

where we have kept up to fourth order in the memory 
kernel and defined 



£ (4) (s) = £< 4) (s + iil n + iAil). 



It is often useful 2 ^ to remove the high-frequency oscil- 
lations arising from the prefactor e H n n+An)i in Eq. 
by introducing a "co-rotating" coherence measure x: 



x(t) = 2e-' l(n - +A ^ t (S + (t)). 



(66) 



Note that £^ 4 '(s) is the fourth-order term of the Laplace 
transform of the memory kernel in the rotating frame: 



£(i) = e- 4 (°"+ A °)*I](t). 



(67) 



This quantity serves as the kernel in the integro- 
diffcrcntial equation governing the evolution of x(t): 



b(t) = -iAilx(t) - i [ dt'Y,(t - t')x(t'). 
Jo 



(68) 



The rotating frame renders some of the more subtle fea- 
tures arising from the hyperfine flip-flops more transpar- 
ent, and we will make extensive use of it when we solve 
the uniform coupling model using the NZ GME later on. 



D. Fourth-order memory kernel 

The fourth-order terms which emerge from inserting 
Eq. ((42]) into Eq. (gBJ are 



(65) 



£( 4 '(s) = iTr{S + P[l-iL Q QG Q (s)}L v QG Q (s)L v Q 
x G Q (s)L v QG Q (s)QL v PS- PI (0)}. (69) 

We can again use Eq. (|4"5j) to replace QGq(s) with QG(s) 
and to discard most of the projection operators, yielding 

£( 4 )(s) = iTr{5+[l -iL Q QG(s)}L v G{s)L v Q 

x G{s)L v G{s)L v S- Pi {0)}. (70) 

If we then replace each instance of the propagator G(s) 
with J dte~ st G(t), then the remaining steps can be per- 
formed in a manner quite similar to the treatment we 
have given for the second-order memory kernel. This 
time, however, the expressions are much more compli- 
cated, and we relegate them to Appendix lAl 



E. High-frequency limit 

From the explicit expression for the fourth-order mem- 
ory kernel given in Appendix [Al it is clear that the dom- 
inant contribution will come from the region s ps — iil n . 
This is because the fourth-order memory kernel is a sum 
of many terms where each term is a product of three 
simple poles. Some of these poles are located at values of 
s which only depend on the fcth nuclear Zeeman energy 
and hf coupling. However, most of the terms contain 
a pole located in the vicinity of s ~ —iil n . Due to this 
structure, x{s) has poles at both s « — iil n and at low fre- 
quencies, but the residues at the latter poles are strongly 
suppressed compared to the former ones. This means 
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that if we consider the function £( 4 )(s) = £( 4 )(s + ift n ), 
we can neglect s, u>k, and A k relative to ft„. The resulting 
approximate expression for E^ 4 '(s) is given in Appendix 
|A"1 The approximate expression for S^ 4 ^(s) can be simpli- 
fied by using the results of Appendix [B] to evaluate bath 
correlators, with the result 



£(4) 



lens E A *4E^ 

f 1 



(t)- (t)+ 2s + ^(^fc ~ ^l) 

Zk ° l s + i{A k -A t ) 



s - i(u)k —ut— h{Ak - At)) 



+ i(u k -ue + h{A k -A e )) 



16ft?s 



i \ k 

5>2 [ c * 




(71) 



The symbols (4 are defined in Appendix [B] The 
fourth-order memory kernel given in Eq. (|7ip simpli- 
fies considerably if we assume a uniformly polarized nu- 
clear bath. In this case, the last two terms in Eq. (|TTj) 
cancel each other. This follows from the fact that 
J2i Pn^k ± ^i >± = c±c± f° r sucn a bath; this identity is 
proven in Appendix |Bj (Here, cj; = is independent of 
k by definition.) In using this form of the identity, we are 
also assuming a homonuclear bath since we have thrown 
away the species information in discarding the site indices 
k and t. This assumption could easily be lifted by intro- 
ducing additional indices, but we will not do this for the 
sake of simplicity. Therefore, for a uniformly polarized 
homonuclear spin bath, we have 



g(4) 



4ft2 



E 



+ i(A k -A e ) 



(72) 



Implementing the same procedure for obtaining the 
high-frequency limit on the second-order memory kernel, 
Eq. ([59]) . we find that the result is a constant: 



V. THE NZ SOLUTION FOR THE UNIFORM 
COUPLING MODEL 

We proceed to solve the uniform coupling model within 
the NZ framework in the case where all nuclei have spin 
1/2. We will find that the solution is structurally incom- 
patible with the exact solution reviewed in Section IIII1 
except at very early times. Let us define the characteris- 
tic time 



n n N 

A A 



(74) 



Note that r Aft" 1 . For times t < r, the NZ solu- 
tion agrees quite well with the exact solution; however, 
we will show that beyond this time scale, the NZ result 
rapidly breaks down. We will argue that this failure is 
an unavoidable consequence of the basic structure of the 
perturbative NZ approach with projector as defined in 
Eq. (|3"Tj) , suggesting that this approach is inappropriate 
in the case at hand. A remedy for this problem will be 
presented in Section IV CI 

Starting from the expressions for the second and 
fourth-order memory kernels, Eqs. ((59)) and (IA21|) . as- 
suming a homonuclear bath, lj^ — w, and setting A k — 
A = A/N, we find in the rotating frame 



IflC 



s + iwi 



(75) 



£(4) 

with 



VUJ\ 



W3 



ITU 2 



2i/i c + c ifi c + c 



ITU?, 



(76) 



n n + An-u-A/(2N), 
n n + An-uj + A/(2N), 

2ft„ - 2lu + Aft, 
Aft, 



and 



A*_ 

4A 



T 



(77) 



(78) 



In the limit of uniform couplings with spin 1/2 nuclei, 
the Lamb shift reduces to 



g(3) 



4ft„ 



E^ 



-Aft. 



(73) 



In the last equality, we have pointed out that the high- 
frequency limit of — is just the Lamb shift we have 
already computed. This was given in Eq. (RkJ)) in the con- 
text of a more general bath (not necessarily uniformly 
polarized). We see that the Lamb shift approximation is 
automatically incorporated into the high-frequency ap- 
proximation. 



Aft 



A 2 



JL 

£2 ) 1 



(79) 



In the above results, we have used the fact that we are 
restricting to the case of a uniformly polarized nuclear 



bath so that the bath correlators 



± 



are indepen- 



dent of the nuclear site index k. Further details about 
how the dependence on arises in (|T5[) and (l76l) arc 
given in Appendix |B] For the remainder of this section, 
we will absorb u> into the definition of ft„. In the rotating 
frame, x(t) must depend on w only in the combination 
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Q — u, so it is simple to restore the explicit w-dependence 
if desired. In any case, it is generally safe to neglect uj 
since for any finite magnetic field w <C f2. 

The solution to Eq. ([68| for the electron spin coherence 
x(t) in terms of a Bromwich inversion integral is 



i r + 



7— too 



dse st : = . 

s + iAQ + jS( 2 )(s)+ffi( 4 )(s) 

( 8 °) 

where xq = x(t — 0) is the initial condition in the time 
domain. Note that we are not making the Lamb shift 
approximation here since we are keeping the full E^ 2 ^(s). 
The Laplace transform of x(t) is a rational function: 



x(s) = x 



*(*)' 



(81) 



where the polynomials R(s) and $(s) are given by 
R(s) = (s + ici7i) 2 (s + iw2) 2 {s + iws)(s + ivja), (82) 



$(s) = (s + iAft)R(a) + fic 



1 + iw± 



jJLC 



Rjs) 

• + IVJ2 



— /.i 2 c + c [2s + i(wi + c^)] 2 [3s + ivus + 2ivo^\ . 

(83) 

The fact that x(s) is a rational function guarantees that 
the sum of the residues equals Xq so that Xq is indeed the 
initial value of x(t). Denoting the seven zeros of &(s) by 
Si , we may write 



$(s)=n( s - s 0- 



(84) 



We can solve for the Sj perturbatively in the small param- 
eter -^p- <C 1 (equivalently, fl n r 1). This essentially 
amounts to an expansion in small fj, about the zeros of 
sR(s), which are 0, — iw\, —IW2, —ivj 3 , —ivo±. Since 
—ivj\ and — iw2 are each zeros of sR(s) with multiplicity 
2, they will each give rise to two separate zeros of $(s). 
We must also keep in mind that the Lamb shift is linear 
in fj, (see Eq. ([75])). We have 



rjfic 



g(f) 



(85) 



2 + - 
H c c 



[2s + i(z&i + T&2)] [3s + 1073 + 2iwi\ 



We have also introduced the parameter 77 so that we may 
consider the effect of making the Lamb shift approxima- 
tion, which would amount to neglecting both the Af2 
and T>( 2 \s) (under the assumption that they cancel one 
another) in the denominator of Eq. (I80|) . This approxi- 
mation is implemented by setting 77 = 0. 



Having introduced all the bookkeeping devices we will 
need, we can proceed to compute the zeros of ^(s) ap- 
proximately by finding the Ui which solve the equation 



^(s^i + fiUi) = 
to leading order in fi, where sr^ satisfies 

SR,iR(sR,i) = 0, 
Following this recipe, the Sj are found to be 
1 r 



(86) 
(87) 



so 



2r 



1 + VI + 16c 



1 - v 7 ! + 16c 



S'2 = —IW\ 



S3 = — IW\ 



It 



2t 



-IT] 



\J Ac + Jc — 77 



-if] + y/ 4c+ /c — T] 



s 4 = -ITD2 + — 
It 



—ii] — yAc /c+ — 77 



S5 



-1W2 



2t 



-irj + y/ Ac /c+ — 77 



s 6 = -im 3 . 



We have only kept up to first order in A 2 /NQ^ = 
(J7„r) _1 in the above expressions for the Sj. Notice that 
sq = —ivDs, so that this root of $(s) is also a root of R(s) 
and thus is not a pole of x(s). The residues of x(s) at 
the poles (s, for i = 0...5) are 2mri where 



Xq I 77 \ 1 — + 16c+C 

r o = — I -1 

c^c 



tt n Tj Vl + 16 



aw, ?7 \ 1 + VT+ 16c+c- 



2 V VnTj VI + 16C+C 



r 2 = -7^0 



r 3 = ix 



r4 = —ixo 



c - ^77 + ^4c+/c - 7/) 
4ft„r ^4c+/c- -77 

c - + y/Ac+/c- - 77^) 

40„r ^4c+/c- -77 

+ ^7/ + \jAc~ /c+ - 77^ 



4ri„r 



nT y/Ac~ /c+ - rj 
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r 5 = ix - 



z + (—ir) + ^/4c /c+ — ^ 



(89) 



40 n r ^4c"/c+ - ?7 
The electron spin coherence in the rotating frame is then 



i=o 



(90) 



A. Time scale for validity of the NZ GME uniform 
coupling model solution 

Notice that two of the exponentials contributing to x(t) 
always diverge in the large-time limit. The terms asso- 
ciated with the poles S3 and S5 exhibit this divergence 
for nearly all nuclear polarizations. The only exception 
is the case of maximal polarization, for which c + c~ = 0; 
we will return to this special case shortly. The fact that 
divergences arise suggests that the NZ result for the uni- 
form coupling model breaks down on a timescale given 
roughly by r ~ NVL/A 2 . The positive real parts of the s, 
are direct contributions of the fourth-order memory ker- 
nel. Even though these contributions appear to introduce 
a time scale cutoff for the solution, their presence actu- 
ally extends the validity of the NZ result to larger times 
as will be shown in the next section when we compare 
the NZ solution with the exact solution. 

Although the large-time divergences generated by the 
positive real parts of the Si indicate that there must be a 
time scale cutoff, they are not really responsible for this 
cutoff. In particular, even if no divergences were present, 
the result we have obtained for x(t) in the uniform cou- 
pling model would still be invalid beyond the time scale 
of t. This is because we are solving perturbatively for 
the phases Im(si) in the phase factors exp[ilm(si)f] ap- 
pearing in Eq. (|90[) . If we solved for Im(sj) to order 
A 2 /NVt n ~ 1/t, then the result would be valid up to 
time scales which are inversely related to the next order 
in the flip-flop expansion (t ~ N 2 ^!^/ A ) since correc- 
tions to the periodic phase factor exp[Tm(si)i] resulting 
from a small correction Ssi will become significant when 
lm(Ssi)t ~ 1 regardless of the fact that A/\fN£l n <C 1. 
However, it is important to note that we have not ob- 
tained the full 0(A 2 /NQ n ) contributions to the Sj as is 
suggested by the particular form of the expressions given 
in Eq. (|88j) . Specifically, if we focus only on the fourth- 
order memory kernel contributions by setting rj = 0, it 
is evident that the Si still receive corrections at order 
A 2 /NQ n ~ 1/t. Therefore, we see that higher-order 
terms in the flip-flop expansion modify the lowest-order 
corrections to the Si, and we would need to include all 
terms in the expansion of the memory kernel in order 
to obtain the full 0(A 2 /NSl n ) corrections to the s». Of 
course, keeping all these terms would yield the exact s» 
and Ti and not just the 0(A 2 /NSl n ) corrections. We 
conclude that the NZ approach with the simple projec- 
tor defined in Eq. (1311) and in which the memory kernel 



is computed up to some finite order in flip-flops will yield 
a solution for the uniform coupling model which is valid 
only up to the time scale of t < r. We have checked that 
even if one computes the poles of $(s) exactly (numeri- 
cally), the NZ result still breaks down at the time scale 
t ~ t, in further support of the conclusion that the break- 
down on this time scale can only be avoided by keeping 
higher-order terms in the memory kernel expansion. 

A notable exception to the above argument arises in 
the case of maximal nuclear spin polarization m = ±JV/2, 
where either c + = or cT = 0. In this case, the fourth- 
order memory kernel vanishes identically and so does not 
modify the Si, suggesting that the NZ result is valid on 
arbitrarily large time scales. We will see below that the 
NZ approach incorporating only the second-order mem- 
ory kernel yields the exact solutionis in this special case, 
confirming this conjecture. 

Note that the above conclusions are compatible with 
the considerations of relevant timescales in a realistic, 
inhomogeneouly coupled system. There we expect the 
box model to apply at short times, t <C N/A. The pre- 
vious works applying the NZ theory to such a problem 
were focusing on the high field regime, A/Q < 1, in which 
the long-time solution could be controlled perturbatively. 
Under this condition, t>N/A, and the NZ theory suc- 
cessfully describes the short-time coherence dynamics. 28 
At low magnetic fields, however, we have t < N/A, and 
the NZ approach is not expected to work properly at the 
short timescale. 



B. Comparison with the exact solution 

When the nuclear spin polarization is maximal, m = 
±7V/2, the fourth-order memory kernel vanishes identi- 
cally. For concreteness, we consider the case m = N/2 
(c + = 0, c~ = 1), in which the Laplace transform of the 
spin coherence becomes quite simple: 



x{s) = x - 



IVJl 



sis + vuj\) 



(91) 



Since there are only two poles, we can obtain these ex- 
actly, 



V3\ % I 9 

s ± = ± -yrof + 4/i, 



(92) 



and it is likewise a simple matter to compute the residues: 

ZD l 



2mr± = — 



1 ± 



The electron spin coherence is then given by 



(93) 
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Plugging in the expressions for w\ and [i, Eqs. (|77p and 
(|78p. we obtain the exact solution from Eq. (|25p with 
to = AT/2. The fact that the second-order memory ker- 
nel suffices for obtaining the exact solution^ for any dis- 
tribution of Ak in the case of maximal polarization was 
previously pointed out in Ref. 0. 

Returning to the case of arbitrary nuclear polarization 
to, it is straightforward to check that the NZ result agrees 
with the exact solution for short times. In particular, if 
we consider the case t -C r, then we may approximate 
the poles as follows 



so 



0, s 2 w s 3 



-IZJi, S4 SS S5 



leading to the short-time behavior 
x 



it) 



x 



fL,T 



(95) 



(96) 



Now turning to the exact solution, Eq. (f2"5"|) . we can 
implement the same approximation by replacing N^ m as 
CC7 2, NJ m ~ w\ in the arguments of the sines and cosines 
arising from a* m and Cj m . Furthermore, we expand the 
coefficients of the sines to second order in the hypcrfinc 
coupling: 



7+ 
fjrn 

N+ 

Z jm 
N im 



1 



A 2 



2zu 2 ~N 2 



A 2 



2vj{N- 



+ 1) — m(m + 1) 



\j{j + l)-m{m-l)}. (97) 



At least for low to moderate nuclear polarizations and 
considering that for typical values of j we may write j(J+ 
1) ~ N, this approximation is valid roughly in the limit 
£l^S>A/ yN. In the above expressions, we can replace w\ 
and W2 with Q n since the differences lead to higher order 
corrections in A/(CIN). After some algebra, we arrive at 



x(t) 



x <l 



N/2 



n nT Nz ^ 3 

j=\m\ 



[jCj + i) 




+ 



x 

n n TNZ" 



N/2 



n ri rNz 



E n j UU + 1) - m ( m + 1)] 

j = \m\ 
N/2 

e -feii n, + 1) - to(to - 1)] 

j=\m\ 



The identity 



2 



N/2 

E 

j = |m 



«j + 1) - ™ 2 ] = i) 



immediately implies that (using also Eq. (|B7I) 

TV/2 

^ E n J D'^' + !) - ™ 2 ± H = cT > 



(99) 



(100) 
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and m = 



1: Exact solution of the uniform coupling model (Eq. 
vs. NZ GME result (Eq. {50)) for A = O, w/fi = 10~ 3 , 
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FIG. 2: Zoomed in version of Fig. [T]along with the Lamb shift 
approximation (LSA) of the NZ GME result for AT = 100. 



and we get back the short-time NZ result, Eq. (|96|) . 

In Fig.[TJ we compare the NZ result with the exact solu- 
tion in the case of zero net nuclear polarization (to = 0), 
and it is clear that the two agree only for short times 
t < t/2. We have not specified the number of nuclei in 
the caption of Fig. [1] because the curves are valid for a 
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FIG. 3: Zoomed in version of Fig. Q]along with the Lamb shift 
approximation (LSA) of the NZ GME result for AT = 10 4 . 
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FIG. 4: Exact solution of the uniform coupling model (Eq. 
(25)) vs. NZ GME result (Eq. O) for N = 10 4 , A = fi, 
cj/fi = 1CT 3 and m = N/4. 
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FIG. 5: Zoomed in version of Fig.[4]along with the Lamb shift 
approximation (LSA) of the NZ GME result for N = 10 4 . 



A/SI between 0.1 and 10, the envelopes in Fig. Q] do not 
change, only the fast oscillation components change; the 
(dis)agreement between the NZ solution and the exact 
result remains unchanged. 

Both the envelope and modulation change under varia- 
tions of the nuclear polarization to. The case m = N/4 is 
shown in Figs. 2] and [5j Fig. 2] reveals that m = N/4 is in 
some sense special because the NZ approach is able to re- 
produce the oscillation period of the x(t) envelope along 
with the amplitude of its real part. We have tried various 
nuclear polarizations and found that the NZ approach is 
only able to reproduce these features in the to = N/4 
case. Generically, the single-projector NZ result visibly 
disagrees with the exact solution beyond t ~ r. 

The fact that Figs. QJ5] show that the NZ solution, 
Eq. (|90"1) . disagrees with the exact answer in Eq. (|25|) 
beyond time t ~ r is not surprising in light of the fact 
that the exact solution is a sum over a large number of 
oscillatory functions, while the NZ result consists of a 
sum of only six exponentials. In the case of a realistic 
wavefunction, the NZ approach can lead to more com- 
plicated behavior due to an appearance of branch cuts 
in the continuum limit&2& It is should also be strongly 
stressed that for a realistic wavefunction, a new pole re- 
lated to flip-flops between nuclei with similar, but dis- 
tinct, values of Ak, appears. In fact, this pole dominates 
the high-field and long-time decay of coherence, and its 
influence is well-controlled for times much longer than r 
and T 2 w r(^), as discussed in Ref. [M This is in stark 
contrast to the behavior of a few poles determining the 
coherence dynamics in the uniform coupling case. 



wide range of N. This is because the envelope of x(t/r) 
is essentially independent of N as can be seen by noticing 
that the envelope arises from the first two contributions, 
roe s °* + rie sit , and these contributions only depend on N 
through their r-dependcnce^ The remaining four con- 
tributions to x(t) give rise to a small modulation which is 
visible in Figs. [2] and [3] The frequency and amplitude of 
this modulation are given roughly by £l n T and l/(fi„r) 
respectively. Since O n r = 4Nfl^/A 2 , we see that the 
frequency scales linearly with N, while the amplitude 
scales like 1/N. This behavior is evident in a comparison 
of Figs. Hand ED 

It should be stressed that while the value of A/ ft plays 
a crucial role for the long-time results in the non-uniform 
coupling case^2^ this parameter does not have any par- 
ticular significance in the uniform coupling case consid- 
ered here. The amplitude of the fast oscillation in the 
exact and NZ results depends of course on the ratio of 
A/Vt, but for large N the amplitude of this oscillation 
does not become significant in the range of which we 
consider, i.e. A/VN. On the other hand, the en- 
velope of the NFID signal depends on fl and A only 
through t ss {N/A)-{Q./A). When varying the parameter 



C. The NZ solution of the uniform coupling model 
using correlated projectors 

We have seen that a straightforward application of the 
standard NZ GME to the uniform coupling model yields a 
result which is in good agreement with the exact solution 
only for very short times. From the technical point of 
view, this failure of the NZ approach could be traced back 
to not having enough poles in the Laplace transformed 
solution, i.e. not having enough frequencies to sum over 
in the real-time expression. Clearly, in order to reproduce 
the exact answer one needs to find a way of bringing all 
these missing frequencies back into the theory. 

The above NZ results were obtained using a standard 
choice for the projection operator P, which was defined 
in Eq. (I5T|) . However, it was shown in Ref. that the 
standard projection operator is far from being the best 
possible choice in contexts where the Hamiltonian ex- 
hibits a significant degree of symmetry. When symme- 
tries are present, one can instead replace P with a series 
of so-called correlated projection (CP) operators which 
project onto invariant subspaces of state space, enabling 
one to expand the reduced density matrix for the system 
as a sum of matrices, each capturing the components of 
the state lying in a particular subspace. This greatly 
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enhances the number of dynamical degrees of freedom, 
resulting in a remarkable improvement in the agreement 
between the NZ calculation and the exact result. In this 
section we will show that, whereas the standard fourth- 
order NZ result agrees rather poorly with the exact solu- 
tion, even at second order the NZ GME agrees extremely 
well with the exact answer when correlated projectors are 
employed. 

For this calculation, we choose to work in the inter- 
action picture defined with respect to the unperturbed 
(Zeeman and Overhauser) part of the Hamiltonian, Hq. 
In terms of the Liouville operator Lj and total density 
matrix p in the interaction picture, the second-order NZ 
GME equation for the electron spin coherence /5 e , |_ be- 
comes 

h,-+ = -l dt'Tr{S + PL I {t)L I (t')Pp{t / )} . (101) 
Jo 

For the uniform coupling model, Lj acts on an arbitrary 
density matrix as 



algebra 



,-iZ+(t-f) ( X +J2 



(106) 



where X^ m and were defined earlier in Eq. ([21 



Taking the Laplace transform of both sides of Eq. (|106l) , 
we obtain an algebraic equation for the Laplace transform 
of y5g™ + (t) (which we call p(s) jm for brevity) which is 
readily solved: 

i? m {s) = pt(s + iZ+)( S -iZ-)L s + i(Z+-Z-) S 2 



(107) 



+ [ZmZ m + -(X+ m ) 2 + -(X jm ) 2 }s 



P++ P+- 
P-+ P- 

1 / h + -p- + - p+-h_ + h+-p— 

2\h-+p ++ -p — h-+ h-+p+- 



- P++K- 

- p-+h+- 



(102) 



with 



= ( A / N yV-A/(2N)) tl - e i{A/N)tI^ 

h-+ = (A/N^-^+^V^I+e-W^ 1 * = h\_, 

(103) 

where I ± = J2k are ^ ne total nuclear creation and 
annihilation operators. At this point, we are ready to 
define the correlated projection operators appropriate for 
the uniform coupling model. In Section Mil we already 
made use of the fact that in the case of uniform couplings, 
the total nuclear angular momentum is conserved, and we 
may work in the \jm) basis of nuclear states. Defining 
the nuclear operator Tlj m to be the projector onto the 
subspace of fixed j and m quantum numbers, we choose 
the superprojector P to be such that it acts on the total 
density matrix as2£ 



pp = ]T Tr z (n jm p) ® — n im = J2 pT ® — 



n 



jm- 



(104) 



The p> e m are a set of matrices which sum to give the 
reduced density matrix for the electron spin: 



N/2 N/2 

Pe=J2^= E E^pr- 



(105) 



jm m= — N/2 j=\m\ 

Inserting Eq. (I104[) into Eq. (jlOll) . we obtain after some 



Pq™ is the intial value of p^L + {t). To obtain fP^L+ii), 
we must use the Bromwich inversion formula: 



~irn 
Pe- + 



py+ioo 

(t) = ^- Jim / dse st ? m ( S ). 



(108) 



The value of this integral is determined by the three 
poles of p jm (s), which we will call sj m , s^" 1 , s J 3 m , and 
by their associated residues r^ m , r^" 1 , rj™. In terms of 
these poles, the solution is 



i=l 



r? m e s - * 



(109) 



where the residues are given by 
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To complete the solution, we must give the explicit forms 
of the si m : 
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FIG. 6: Exact solution of the uniform coupling model (Eq. 
(25)) vs. NZ-CP result (Eq. JTTUl) for TV = 10 4 , A = ft, 
w/ft = 10" 3 and m = 0. 



FIG. 7: Zoomed in version of Fig. [6] 
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Since the interaction picture can be thought of as the 
frame rotating with frequency O n , to switch to the ro- 
tating frame with frequency O n + AO we simply need to 
multiply the result by e _iAa *. The electron spin coher- 
ence in the rotating frame and in the presence of fixed 
nuclear polarization m is then 



XNZ-cp(t) = -g-e 
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(114) 



This solution is plotted along with the exact solution in 
Figs. [5] and [7J It is clear from these figures that the 
NZ solution with correlated projectors (NZ-CP) agrees 
remarkably well with the exact solution. 

It is not difficult to extract analytically the envelope 
and modulation from the NZ-CP solution. To do this, it 
helps to first perform a large N expansion, wherein the 
leading order terms in the s^ m are 
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In the limit of large N, the coherence therefore takes the 
form 
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(117) 



The factor in curly brackets in Eq. (|117[) is precisely that 
which gives rise to the small modulation depicted in Fig. 
[Jj while the rest of the expression produces the envelope 
in Fig. |U 

The above expression can be obtained from the exact 
solution, Eq. (|25|) . by expanding in the limit X^ m Z m , 
in which case we find 
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where a jm , Cj m and TVj^ were defined in Eqs. ([20)1 and 
(H^) ([2T1) . These expressions result from expanding only the 
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factors appearing in a jm and Cj m . If we further- 
more expand the N- m factors in the temporal exponents 
as 



(x 



A\Z?, 



(119) 



and also make the following approximation in the denom- 
inators, 



7 d 



(120) 



then we get back Eq. ([TT71) . 

If we assume typical values of j and m with j > m, 
then we have roughly X^ m ~ A/ VN, and the above ap- 
proximations are valid in the limit Q ^> A/VN. Recall 
that this is the same condition we had for the validity of 
the NZ solution without correlated projectors. In fact, 
the approximations we have made in the exact solution 
to arrive at the NZ-CP result are essentially the same 
ones we made to relate the exact and NZ solutions (see 
Eq. ([97]) ). the only difference being that here we have 
kept some of the j-dependence in the exponents when 
evaluating the sum over j. Moreover, note that the ex- 
pansion of the Nj m factor appearing in the temporal ex- 
ponent of the exact solution introduces a time scale on 
which the solution is valid. It is somewhat difficult to 
say precisely what the time scale of validity is, however, 
because of the sum over j. To obtain a very rough esti- 
mate, we can again consider a typical value of j and write 
Xf m ~ A/y/N. The timescale should roughly correspond 

to the inverse of the error introduced by expanding Nj^ 
in the temporal exponent. This error is on the order of 
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(assuming that j 3> m), indicating a time scale 
t ■ 
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D. Comparison with the RDT solution of the 
uniform coupling model 

We now review the effective Hamiltonian RDT solu- 
tion to the box model and show how it is related to the 
exact solution and the NZ-CP result of the previous sec- 
tion. At second order in the flip-flop expansion, we have 
processes in which the electron spin flip-flops with the 
kth nucleus and then subsequently flip-flops with the £th 
nucleus. This leads to an effective flip-flop interaction 
between the fcth and ^th nuclear spins which is medi- 
ated by the hyperfine interaction. We can derive the ef- 
fective Hamiltonian describing this internuclear flip-flop 
process by performing a canonical transformation on the 
hyperfine Hamiltonian in Eqs. ©-([S]): H e ff = e~ s He s , 
where the unitary operator e~ s is chosen in such a way 




FIG. 8: Exact solution of the uniform coupling model (Eq. 
HH)1) vs. RDT result (Eq. (TTJU) for TV = 10 4 , A = Q, 
Lj/Q = 10~ 3 and m = 0. 



as to eliminate the original Vg interaction ! 11 ' 12 ! 23 ' 26 It 
is important to note that the transformation on states, 
— > e~ s \ip) 1 is neglected in the derivation of H e fj. 
The dominant terms in the resulting effective Hamilto- 
nian are 
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The effective Hamiltonian solution for NFID was ob- 
tained in Ref. [26| in the case of zero nuclear polarization, 
m = 0. In the rotating frame, it is given by 
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(124) 



This result is compared with the exact solution, Eq. (|25[) . 
in Figs. [S] and [S] To generate these figures, we have 
used the high-frequency form of the Lamb shift, Eq. (|63l) . 
where for zero polarization, = 1/2 (see Appendix IB1 . 
It is clear from Fig. [8] that the envelope of the RDT so- 
lution agrees very well with that of the exact solution up 
to time scales well beyond r. The fact that the two en- 
velopes agree even at very large time scales is particularly 
interesting given that, in the case of RDT, x ~ 1/t for 
t t; this 1/t tail is evidently reproduced by the sum 
of many exponential terms comprising the exact solution. 
On the other hand, Fig.|9]reveals that RDT does not suc- 
ceed in capturing the fast small-amplitude modulation of 
the coherence. 

It is not difficult to extract the RDT solution as a limit 
of the exact solution. First recall that the RDT result 
using the effective Hamiltonian is expected to be valid 26 
for Vt ^> A/y/N and in the limit of large N. Starting 
from the exact solution, Eq. ([25]) . we can apply the first 
of these approximations along with the fact that for the 
j states which contribute the most to the sums we have 
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FIG. 9: Zoomed in version of Fig. |8] 



factor appearing in the temporal exponent of the exact 
solution that introduces a time scale on which the NZ- 
CP solution is valid. The fact that we are making the 
very same expansion here implies that the RDT solution 
will be valid on a time scale similar to that of the NZ-CP 
result. 



VI. DISCUSSION 

In this Section, we aim to understand the reasons 
for the failure of the standard (single-projector) NZ ap- 
proach discussed in Sec. |V] and for the success of the 
correlated projector approach from the previous Section. 



j(j + 1) ~ N, to write 



7 ± 



N 



A 1 



2n 2 N 2 



j(j + 1) 



so that we obtain 



N/2 



, (125) 



(126) 



Expanding the exponential and using that in the large N 
limit, 
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we arrive at the RDT solution 
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Strictly speaking, this derivation is only valid for t < r 
since otherwise the infinite series in Eq. (|128|) is not de- 
fined, and the last line in Eq. (|128|) should be treated 
as an analytical continuation of the previous expression 
along the real axis (a similar kind of analytical contin- 
uation was in fact encountered in the derivation of the 
RDT result in Ref. 261). The correctness of the analyt- 
ical continuation can be easily checked by comparing a 
numerical evaluation of Eq. (|126|) with Eq. (|124l) in the 
limit N — > oo for all times t. 

Notice that the approximations we have made in the 
exact solution to arrive at the RDT result are essen- 
tially the same ones we made to relate the exact and 
NZ-CP solutions (see Eq. (|119[0 . the only difference be- 
ing that in the context of NZ-CP, we kept one more order 
in the expansion, yielding the additional factor in curly 
brackets in Eq. (|117[) . This additional factor is precisely 
that which gives rise to the small modulation depicted in 
Fig. Thus, the RDT and NZ-CP solutions of the box 
model belong to the same expansion of the exact solu- 
tion. Moreover, note that it is the expansion of the Np m 



A. Non-narrowed Gaussian bath 

It is instructive to take a small detour, and discuss first 
the NZ solution for the case of Ising coupling to the bath 
spinsj 5 -^ for which Vg = 0. Furthermore, let us abandon 
for a moment the assumption that the state of the bath 
is narrowed, and take pi ~ 1. 

This case is of course exactly solvable^ and in the 
large N limit the results are equivalent to the classi- 
cal calculation,— in which the bath is assumed to be 
static, and the central spin dynamics is calculated by 
averaging the spin precession over a Gaussian ensemble 
of nuclear polarization vectors Bj of the form P(B|) ~ 
exp (-B 2 /2a 2 ), where a 2 = ± £ fc I k (I k + l)A% A sim- 
ple Gaussian integral gives the well-known exp[— (t/T 2 *) 2 ] 
decay of the envelope of the FID signal, with T 2 * = y/2/a. 
For large N, this is practically equivalent to the evalua- 
tion of the quantum expression 



<S+(i)> = ^e^^exp( ? ^i)(5 + (0)) 



(129) 



Note that in the above expression, every \m) state present 
in p/(0) contributes a frequency h^., and the averaging 
of these frequencies leads to dephasing (inhomogeneous 
broadening) . 

It was noticed before that the standard NZ approach, 
when carried out to a finite order in the expansion (in 
powers of h z in this case), cannot reproduce the above 
simple solution. 59 On the other hand, the TCL approach 
gives the exact result already in the 2nd order of expan- 
sion. This feature can be traced back to the fact that 
in the NZ expansion, one obtains so-called partial cumu- 
lant o 43 ^ 59 of the bath operators (i.e. (h^) pc = (hf) — (hi) 2 
appears in the 4th order), while in the TCL calcula- 
tion one encounters the ordered cumulants (i.e. (h A z ) oc — 
(h A z ) — 3 (hi) 2 ), which are clearly closely related to the 
usual cumulants. For the Gaussian bath, the ordered cu- 
mulants beyond 2nd order are identically zero, and the 
TCL approach succeeds because it captures the crucial 
statistical properties of the bath already in the 2nd or- 
der of the expansion. 61 The NZ approach, on the other 
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hand, is incompatible with the structure of the bath cor- 
relators, and it has to be carried out to infinite order to 
recover the Gaussian decay of the transverse spin. 



B. Uniform coupling model as a classical 
non-Gaussian bath 

The effective Hamiltonian-based RDT solution for de- 
phasing in the uniform coupling model, Eq. (|124[) , which 
is the same as the large N limit of the exact solution (see 
the discussion leading to Eq. (|128p ). can also be obtained 
from a classical calculation.— We simply fix the B\ com- 
ponent of the nuclear field, perform the Gaussian average 
of the classical equations of motion for S x,v , and recover 
the RDT results. The fast oscillation of the exact solu- 
tion is clearly an effect which is being missed by the clas- 
sical large- N limit. The same result is obtained when one 
employs the classical limit of the effective Hamiltonian, 
S z \(Bf) 2 + (B v I ) 2 ]/2Vt (as was done for the spin echo case 
in Ref. 31), and performs the Gaussian average over all 
the precession frequencies due to the possible magnitudes 
of the transverse nuclear field. This shows that the RDT 
(at short-times in the general case) and the uniform cou- 
pling calculations are essentially equivalent to classical 
averaging over a square of a Gaussian-distributed classi- 
cal variable. The latter approach was succesfully applied 
to the calculation of the Rabi oscillations decay of spin 
qubitsiS£3& 

In the case of a bath operator coupling which is 
the square of a Gaussian-distributed variable, the bath 
correlators are such that cumulants of all orders are 
nonzero. However, these cumulants are simply the ring 
diagrams , 26 ! 44 and the summation of all of them can be 
performed. In order to achieve this, it is crucial to use a 
theoretical approach in which the standard, i.e. ordered, 
cumulants appear in a natural fashion, as it occurs in 
the effective Hamiltonian-based RDT calculation, or in 
the TCL approach^ 3 . 

As we discussed in the previous section, the standard 
(single-projector) NZ theory does not have a natural re- 
lation with the statistical properties of the bath, specifi- 
cally with the structure of the ordered cumulants of the 
bath variables. Furthermore, in the standard NZ ap- 
proach we project on a single tensor product state of the 
system and a fixed state of the environment (chosen here, 
as it is usually done, to be pi(0)). The Ising coupling case 
clearly shows that this can a sub-optimal approximation: 
in this case, the central spin coherence decays only due 
to the fact that the exact density matrix at finite times 
possesses a nontrivial structure deriving from the evolu- 
tion of each of the nuclear states present in pi(0), i.e. in 
Eq. (|129|) we have a sum over all the \rii) states, and each 
state contributes a different frequency. What is explic- 
itly shown in the exact solution of the uniform coupling 
model is the fact that, even for the narrowed state of the 
bath, this is the case. The exact solution in Eq. ([251) in- 
volves a sum over the j quantum numbers: in order to 



obtain the correct result, one has to preserve the struc- 
ture of the nuclear density matrix which is singled out 
by its coupling to the central spin (i.e. the nontrivial dy- 
namics in the j, m, m± 1 subspaces). In the classical limit 
this averaging over the nuclear states present in the nar- 
rowed density matrix is expressed by integration over the 
Bj and B\ components of the nuclear field. The RDT 
calculation captures correctly this averaging by a proper 
resummation of all the diagrams relevant in the large N 
limit. The NZ-CP calculation simply preserves the key 
structure of the total density matrix during the system's 
evolution, and in this way recovers the classical averaging 
limit. 



C. Consequences for the non-uniform couplings 
theory 

As noted in Ref. l22|, the single-projector approach is 
expected to be valid when the system-environment cou- 
pling is weak in some sense. As shown in Ref. 0, and 
|28| . in the case of the hf-coupled spin bath this "weak cou- 
pling" condition is A/Vt< 1, and it apparently guarantees 
that the single projector approach is valid at all times. 
Under this condition the short time (t <C N/A) decay is 
very small, i.e. a quadratic initial decay, x{t) 1 — 2(t/r) 2 
with <<t. This result does not depend on the shape of 
the electron wavefuntion. The decay at longer times is 
due to the fluctuations of the Overhauser field induced 
by flip-flops involving nuclei with different couplings to 
the central spin; the inhomogeneity of couplings is cru- 
cial at this timescale. Most of the decay for .4/57 <c 1 is 
of the exponential form recovered by the RDT calcula- 
tion. With decreasing J7, the non-exponential features of 
the single-projector NZ solution become more prominent, 
but the form of the NFID at low fields f2 < A cannot be 
determined in a controlled fashion, 28 at least for realis- 
tic values of nuclear spin I > 1/2 and for zero or small 
nuclear polarization. 

At short times, the RDT solution remains controlled 
(as long as one believes in the robustness of the trans- 
formation leading to the effective Hamiltonian on this 
timescale) when Q 3> A/y/N. In this regime, the RDT 
solution is explicitly independent of the shape of the elec- 
tron wavefunction ) 25 i 26 as was the case for the short-time 
NZ result i^ 8 . This agrees with the intuition that at this 
timescale all the inter-nuclear flip-flops should be virtual, 
and all of them should be equally important. It suggests 
that on short timescales, it should be possible to map 
the dynamics of the real system onto the dynamics of a 
model system with uniform couplings! 26 ' 36 As we have 
shown here, the RDT reproduces the envelope of the ex- 
act result in this case, while the standard NZ approach 
fails beyond very short times, t <T = 4Nfl/A 2 (at which 
the envelope is well described by the quadratic decay), 
suggesting that the latter approach is indeed unable to 
describe the low-field decay, most of which occurs at short 
times (t < N/A). This is consistent with the results of 
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Refs. I9 U23U281 in which only high fields were considered 
(or a large nuclear polarization was assumed), and the 
single-projector NZ theory was shown to be very hard to 
control at low fields. It should be however noted that 
the broadening of the energy bandwidth available for 
multiple inter-nuclear flip flops in higher orders in Vs, 
discussed briefly in Ref. 28, shows that the boundary of 
the short-time regime can move to even shorter times as 
one goes to a higher order in the expansion of the NZ 
memory kernel. It is nevertheless unclear what happens 
to this trend as the expansion is continued to higher or- 
ders. The previous discussion suggests that in order for 
a single-projector NZ approach to recover the uniform 
coupling bath limit one has to go to an infinite order 
in the memory kernel expansion. This makes it hard to 
perturbatively delineate, within the single-projector NZ 
approach, the timescale on which the uniform coupling 
model gives a good approximation to the real problem. 
A full-Hamiltonian approach starting from an assump- 
tion that the box model holds on some timescale, and 
correctly describing the decoherence in a broad range of 
magnetic fields, should shed more light on these issues. 

The use of correlated projection operator a 22 ' 24 ' 43 
makes the 2nd order NZ result agree very well with both 
the envelope and the small oscillations of the exact re- 
sult for NFID decay in the uniform coupling case. This 
strongly suggests that the correlated projection operator 
technique should be an important element of a theory of 
central spin decoherence encompassing all the regimes of 
magnetic fields and timescales. 



VII. CONCLUSIONS 

In this paper, we have applied the Nakajima-Zwanzig 
generalized Master equation, originally developed for the 
hyperfine-coupled central spin problem in Refs. I9 H23H2H 
to the exactly solvable case of uniform hyperfine cou- 
plings. We have shown that the NZ calculation of nar- 
rowed state free induction decay (NFID), i.e. the evolu- 
tion of the electron coupled to a nuclear bath acting on 
it with a well defined Overhauser field, fails very quickly 
in this case, and can only account for the initial decoher- 
ence dynamics. We have traced the origin of this failure 
to the fact that, in this NZ approach, a single operator 
was used to project the total density matrix onto the 
tensor product of the central spin and bath density ma- 
trices. While this approach works well at high magnetic 
fields, at which the coherence decays at long times, and 
the inhomogeneity of the electron's wavefunction (i.e. the 
inhomogeneity of the hf couplings) is crucial for the co- 
herence decay) 23 ' 28 it fails at short times and low fields, at 
which the shape of the wavefunction should be irrelevant, 
and the uniform coupling model is expected to be appli- 
cable. In the latter situation, one has to modify the NZ 
approach by introducing a family of correlated projection 
operators , 22 i 24 with which one can capture the essential 
features of the exact dynamics: the electron loses its co- 



herence by a rather simple dephasing process, in which 
the electron spin states acquire a different phase when 
interacting with nuclear states from different subspaces, 
and the coherence is lost due to averaging over the ini- 
tial nuclear states. With this modification, the NZ GME 
reproduces very well the exact solution for the uniform 
coupling model. 

The exact result for NFID in the uniform coupling 
model is also reproduced by the effective Hamiltonian- 
based solution from Refs. 25.26 on a timescale long 
enough to capture the full coherence decay. This 
is another example of the effective Hamiltonian-based 
theor y 25 ' 26 reproducing the semiclassical limit (i.e. the 
quantum calculation being equivalent to an average over 
classical nuclear fields) of the central spin dephasing 
problem^ 

According to the effective Hamiltonian-based theory, 
NFID should be described by Eq. f|124[) up to a timescale 
of ~ 10 yus in GaAs dots (with N ps 10 6 ), and this de- 
cay will be significant for magnetic fields smaller than 1 
T. The NFID measurements in this parameter range are 
within reach of the fast measurement techniques devel- 
oped for gated GaAs dots^ Upon increasing the mag- 
netic field, the decoherence times calculated from the 
purely hf central spin Hamiltonian quickly become com- 
parable to the predicted decay times due to dipolarly- 
induced spectral diffusion ; 12 i 52 which were obtained with 
the cluster methods succesfully describing the high-field 
decay of the spin echo signal in GaAs4£ This makes the 
long-time regime (the exponential decay) probably very 
hard to observe in large GaAs dots. On the other hand, in 
the InGaAs dots, especially the smallest ones (AwlO 4 ), 
the short-time regime extends to at most ^100 ns, and 
the low-field requirement might be incompatible with the 
spin splitting needed for the optical manipulation of the 
qubit. These dots seem to be more suited for experimen- 
tally investigating the high-field and long-time regime 
of exponential decay since the predicted 5 -^ decoherence 
times due to the nuclear dipolar interactions are on the 
order of 5-10 /xs (which is a lower bound since quadrupo- 
lar interactions were not included in that theory) , giving 
a window of timescales in which the decay due to the hf 
interaction only could be observed. The exponential de- 
cay was indeed seen in such dotsj 3 -^ but more research is 
needed to ascertain its magnetic field dependence. If the 
nuclear state narrowing could be done in InGaAs dots at 
lower magnetic fields (e.g. around 1 T and possibly be- 
low), the crossover between the short and long timescales 
should be seen, making it possible to check the timescale 
at which the uniform coupling model predictions hold and 
to verify the predictions of Ref. HH in the f2 ss A regime, in 
which the single-projector NZ results differ significantly 
from the effective Hamiltonian results. 

The results of this paper make a strong statement re- 
garding the structure of a theory which could describe the 
electron spin decoherence without the use of an effective 
Hamiltonian in a broad range of magnetic fields. In order 
to properly describe the short time (t^N/A) regime, in 
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which the coherence is expected to decay at low mag- 
netic fields, one should employ the correlated projection 
operator approach in the derivation of a generalized Mas- 
ter equation. This was done employing both the NZ and 
the time-convolutionless (TCL) methods in the uniform 22 
and non-uniform 2 ^ coupling cases (in the latter case only 
using the projectors of spaces of fixed to, not the jm 
spaces used here), albeit only for a thermal nuclear bath. 
In this work we have shown that this feature of the theory 
is crucial also in the case of NFID and uniform couplings. 
The case of low- field NFID in a realistic, inhomogenously- 
coupled system remains to be further investigated. While 
we have focused here on the NZ approach, which was the 
subject of intense researc h^ 23 ' 28 in the context of NFID, 
it might turn out that the TCL approach (untested yet 
for NFID) will be both easier to implement and more nat- 
ural. The latter statement is supported by the fact that 
the structure of TCL is closer to the cumulant expan- 
sion structure which underpins the effective Hamiltonian 
solution from Refs. [25ll26l These theories not only prac- 
tically agree with the NZ approach at high fields in the 
non-uniform coupling case, but they also capture cor- 
rectly the exact solution in the uniform coupling case, 
strongly suggesting that they correctly describe the low- 
field and short-time decay of NFID, just as they correctly 
described the spin echo signaL 31 i 48 However, the nature 
of the crossover between the high-field/long-time, and 
low-field/short-timc decay behaviors remains to be elu- 
cidated using a theory capable of treating both regimes 
on equal footing. 
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Appendix A: The full expression for the 
fourth-order memory kernel 



In this section, we will work out the explicit form for 
the Laplace transform of the fourth-order memory kernel. 
We will not make any assumptions about the distribution 
of hyperfine couplings or about the magnitude of rela- 
tive to other scales in the problem. The only assumption 
is that the initial nuclear density matrix has the form 
shown in Eq. (TT2"1) . Starting from Eq. ((TO)) , we replace 



each occurrence of G(s) with J °° dte st G(t) to obtain 

roc 4 

E (4) (s) = i / T[dt i e- s '£* ti Tr{S + [s-iLoQG(ti)} 
J o t=i 

x L v G(t 2 )L v QG(t 3 )L v G(t 4 )L v S- pi (0)}. 

(Al) 

This step allows us to work with G(t) as opposed to its 
Laplace transform, making it easier to find explicit ex- 
pressions. We will evaluate the string of operators ap- 
pearing inside the trace in two stages. First consider 

T = G(t 3 )L v G(U)L v S- pi (p). (A2) 

We can evaluate this using Eqs. (|55|) and (|56p. finding 



T 



T n 



(A3) 



with 



kf 



lit = ^h^t^Pi^+cuP! {OK It], ( A4 ) 

kl 



where 

a>ki 



e i(u k --£-)t 4 _|_ e i(w f + -/)t 4 
I _ 1 A A e i(Q„+0Jh-ut-^(A k -At))t3 e i(uk+^-)U 
g M — _A £ e i(ttn+0Jk-ut+ji(Ak—Ae))t3 e i(u k -^£-)U 

(A5) 

For the second stage of the evaluation, we define 

R = L v G(t 2 )L v QT, (A6) 



Given the structure of T, Eq. (|A3|) . it is not difficult to 
show that R has a similar form: 



R 



R n 



(A7) 



In terms of R, we may write the trace in Eq. (|A1[) as 
Tr{S + [a-iL QG(t 1 )]R}=Tr{A l Rx}, (A8) 

with 

Ai = s + i(h z - h^ n+h ^ tl . (A9) 

Therefore, it is only necessary to compute one of the 
components of R, and this can be expressed in terms of 
T according to 

Rit = \h+U + (t 2 ){h-T n -h- pi (0)TiiT lt } 

- T n h + }ul(t 2 ) - hj-{t 2 )[h + T n 

- T n h- + Pl {0)Tr[Ttf]h-}uUh)h + , (A10) 
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where U±(t) were defined in Eq. (|54|) . The explicit form 
of the trace, Eq. (|A8[) . is rather messy, so we first break 
it up into four parts in an effort to improve readability: 



Tr{A lJ R lt } = ^X i) 



(All) 



1=1 



x 4 



~Tr {A 1 / l +[7 + (t 2 )%/i+f7|(t 2 )} 

^ , &{A 1 tr_(t 2 )ft+r n ul(t 2 )fe+} , 

-^Tr[T tt ]Tr{A 1 / l +C/ + (t 2 )/ l - /97 (0)C/|(t 2 )} 
jTr[T 4t ]Tr{A 1 f/_(t 2 )p J (0)^-[/i(t 2 )ft + } ) 
jTrjAi/i+tZ+^/i-r^t/l^)}, 
jTr{A 1 C/_(i 2 )2V l -[/i(t 2 )/ l +} . (A12) 



These evaluate to 



-Xi — —j ^ (ikiA p A q 



kipq 



x [s + i(A p - ^) e »(f^+-4p-^)*i 
x Tr{J+7,- pj (0)J^/+}, 



(A13) 



X 2 = e 
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- c k e { ~ t4 + c+e i ~ u 



(A14) 



Performing the four-fold Laplace transform in Eq. (|A1[) 
on each of the Xi and denoting the results by Y i: we find 



Yi = — ^2 A k A e A p A q 



klpq 



Tr {I+Iz P7(0)/-/+} 



s + i(Q n -u k -ue + \ {A k - A f )) 

Ar, - Ag 
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s - i(Q n + A p - At) 
1 



1 



s-i(w p + 4f) s-i(uj q -^-) 



(A17) 
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(A18) 



y 3 = A kAeA p A q 



klpq 
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s - i(ujp + -^) [ 

Tr{/ £ + 47p/(0)/+/-} 



s - i(oj k + 4f) s - i(Q n + uj k - uji - \ {A k - At,)) 

1 Tr{/+/- p/ (0)4~/+} 

s - i{cu k - 4f ) s ~ + Wfc - W£ + \{A k - A t )) 



l + i 



Ap — A q 



s + Ap Aq) 



(A19) 



X 3 = -Y,ApAqe^ + ^\ S b ke Tr{l+I~ Pl (0)I+l;} 



klpq 



+ cm U + i{A p - A q )t 



i(a„+A p -A,j)ii 



Tl-{l+IqPj(0)I^I+} 



(A15) 



Yi = ^^A k A t A p A 
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Tr{l-I+ Pl (0)I^I+} 
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1 Tr {I+I k - Pl (0)lq /+} 
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klpq 



l + i 
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(A20) 



bki U + i{A e - A k )( 



.iidn+Ae-A^t! 



The Laplace transform of the fourth-order memory kernel 
is given by 



Tr{/+7 fc -p 7 (0)J-7+} . 



(A16) 



(A21) 



i=i 
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In the high-frequency limit, we have in the rotating 
frame defined by Vl n (Yl = Yi{s + i£l n )) 



(A22) 



44 + 4] [ej + c+] i, (A23) 
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To^E^ 4 ^ 

Tr{/+/,7p/(0)/+^} 
s - i{uj k -wt- \{A k - 

s Tr{/+/- p/ (0)/,7/+} 

- i(Ap - A q ) s- i{uj k - ut + \ {A k - At)) 



(A24) 
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and 



S( 4 )(s + iQ n ) = E(«) =i^2Yi(a). 
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Appendix B: Nuclear bath correlators 

In this section, we evaluate the nuclear bath correlators 
which arise in the expression for the fourth-order memory 
kernel. For example, in the expression given for Y3 in 
Appendix ([A"]), there appears the correlator 



Tt{I+I- Pi (0)I^I+} 



(Bl) 



Since pi(0) is assumed to be diagonal in the |rxj) basis, 
this correlator is only non-zero when each pair of raising 
and lowering operators act on the same nucleus, and we 
obtain 

Tr{l+I- Pl (0)l-I+} 
= [S k( 5 pq + 6 kp 6 eq (l - S kq )] Pttcg+d®-, (B2) 



where we have defined 



I k (I k + l)-(miy T ml- (B3) 



Recall that I k is the total spin of the fcth nucleus and 
that m\ is the eigenvalue of I k associated with the state 
\rii). (The total spin of the fcth nucleus depends on k 
since we have incorporated information about different 
species into the nuclear site index k.) These quantities 
are related to the correlators according to 



«NTr{/^(0)}=^p4' 



(»)± 



(B4) 



We may express the other four correlators appearing in 
the Yi in terms of the c£ : 

Tr {l+Ilpj(0)I-I+} 

= [SkqSip + S kp Se q (l - Ski)] ^ PiiCk Ci , 

i 

Tr{//4>/(0)/-/+} 
= [heSpg + S kp di q (l - 6 kq )] ^2 Piic£'~ c q ^ + , 

i 

= ]T p u [SuS pq c^-c^- + S kp S eq (l - 5 kq )c®-<f+], 

i 

Tr{l q I+ Pl (0)I-I+} 
= £ PiiW Pq cf + (f + + 6 kp 6 eq (l - 4 9 )4 l)+ c«-]. 

i 

(B5) 

As an aside, we note that in the case of a uniformly 
polarized homonuclear bath, the two-operator correlators 
become independent of k: 



= c t = E p» c . 



k 



(B6) 



For the case of spin 1/2 nuclei, these are particularly 
simple, 



1 \ - , 1 m 
- T ^ Pu m k = - T -, 



(B7) 



where m is the net polarization of all the nuclei defined 
by 



E p «E m *- 



(B8) 



For a uniformly polarized nuclear bath, m is closely re- 
lated to the quantity h z n defined in Eq. (f34|) : 

Nh* 



A 



(B9) 



Next, we will prove a relation which is useful for sim- 
plifying the correlators in Eq. (|B5j) . In particular, we 
want to show that in the case of a uniformly polarized 
nuclear spin bath, we have 
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PiiC k C t — 



E 



Pii°k 



E 



w- 



c c 
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In the final expression, we have used that the condition 
of uniform polarization implies that the result is inde- 
pendent of the nuclear site indices k and £. We have also 
implicitly assumed a homonuclear bath since informa- 
tion about different nuclear species is incorporated into 
the site index. This last assumption can easily be relaxed 
by retaining this information in the form of an additional 
index; however, in the situations in this paper where Eq. 
(|B10[) is employed, the homonuclear assumption is made 
anyway for the sake of simplicity. Also note that we are 
not assuming spin 1/2 nuclei in Eq. ([BIO]) . Replac mg 

the cjj. with explicit expressions from Eq. (|B3j) and 
canceling several terms, we obtain 
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Each term in Eq. (|F311|) is independent of k and I in the 
uniformly polarized case, so that the second and third 
lines on the right-hand side of the equation vanish iden- 
tically. Consider the term J2 i pum k m\ in the last line. 
We can multiply this by A k Ai and sum over k and £. 
Since the \rii) states all have the same h z eigenvalue, we 
know that 



A k A e mimi = {h z n f 



(B12) 



is independent of i. On the other hand, J^i Pii^k m \ * s 
independent of k and £, so that multiplying by A k Ai and 
summing over k and I simply amounts to multiplying this 
expression by (^Z k Akf = A 2 . These two observations 
together imply that 



£ 



From Eq. (|B9I) we have 



E 



P*i m k = T7 = 



m 
N 



K 

A" 



(B13) 



(B14) 



showing that the last two terms on the right-hand side 
of Eq. (|B11|) cancel each other. 



To finish the proof of Eq. (IB10[) . it remains to show 
that the first line on the right-hand side of Eq. (|B1 1|) 
vanishes. This follows straight-forwardly from a simple 
generalization of the previous argument. We begin with 
the following expression 



E 



pum k m e m p m q 



£< 




We can evaluate explicitly each of the two terms by mul- 
tiplying by A k AiA p A q / A and summing over k,£,p,q. 
Since each term is independent of k,£,p,q under the as- 
sumption of uniform polarization, this procedure should 
leave the terms unchanged. We find that the two terms 
evaluate to the same expression and thus cancel: 



E 



™* ™^ ™* 



\ ^ 



pum k m t m p m q - y ^piim k m 




p ii m l m q 



'fi)-(S )=1) . 



The first line on the right-hand side of Eq. (|B1 1|) is just 
a special case of Eq. (|B15I) where p = k and q — £, so 
we have also shown that this line vanishes, completing 
the proof of Eq. (|B10|) . Repeated application of these 
arguments can be used to show more generally that 



E 



{i)± (i)± _ ± ± 



(B17) 



where here the two sets of ± signs are independent of 
each other. 



The identity in Eq. (|B17|) is used extensively in the 
context of the high-frequency limit considered in section 
IIV El We can also use it to simplify the four-operator 
correlators in Eq. (|B5|) . In the uniform coupling model, 
the indices of these correlators are summed over, so we 
have 

E Tr {# J i>(o)J fc - = E Tr {#VM0)4^ + } 

k£pq 



Tr {//7 fe -p7(0)/-J p + } = 2N 2 c + c 



ktpq 



^Tr{/-J+ P7 (0)/ fe -/+}=iV 

k£pq 



2 c+. 



(B18) 



These relations were used to obtain Eqs. (|75t and (1751) . 
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